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Abstract 

In  this  study  we  examine  a  power  spectral  density  metric  for 
characterizing  laser  speckle  and  predicting  its  effect  on  target 
detection.  To  evaluate  the  metric,  we  measured  contrast  energy 
thresholds  on  both  a  laser  speckle  background  and  backgrounds 
consisting  of  randomly  modulated  pixel  luminance  (i.e.,  pixel 
noise).  We  found  that,  at  the  same  power  spectral  density  levels, 
energy  thresholds  obtained  for  gratings  superimposed  on  pixel 
noise  were  slightly  higher  than  analogous  thresholds  obtained 
with  laser-speckle  noise.  In  the  discussion  section,  we  outline 
necessary  improvements  to  the  laser  speckle  power  spectral 
density  metric. 

1.  Introduction 

Simulation  of  the  out-the-window  view  of  pilots  requires  a  high- 
resolution,  wide  field-of-view  display,  whieh  in  turn  requires  a 
high  pixel  eount.  One  approaeh  to  this  problem  has  been  the 
development  of  laser  projeetors  eapable  of  displaying  up  to  20 
megapixels.  One  major  disadvantage  of  laser  projeetors,  however, 
is  the  laser  speekle  that  is  generated  whenever  eoherent  radiation 
is  seattered  from  a  surfaee  whose  roughness  is  eomparable  to  or 
greater  than  the  wavelength  of  the  radiation.  Laser  speekle 
appears  as  a  grainy  pattern  superimposed  on  the  displayed  image, 
whieh  ean  potentially  limit  target  deteetion. 

The  goal  of  this  researeh  is  to  develop  a  metrie  that  eharaeterizes 
laser  speekle  and  quantifies  its  effeet  on  target  deteetion.  One 
metrie  often  used  to  eharaeterize  speekle  is  RMS  eontrast,  or 
speekle  eontrast.  This  metrie  has  the  benefit  of  being  easy  to 
eompute  from  the  speekle  image  in  that  it  is  equal  to  the  standard 
deviation  of  the  image  pixel  luminanee  values  divided  by  the 
mean  luminanee.  However,  speekle  eontrast  is  an  ineomplete 
deseription  in  that  it  does  not  aeeount  for  the  spatial  seale  of  the 
noise  elements  (or  spatial  frequeney  bandwidth)  of  the  speekle 
noise.  A  more  appropriate  metrie  is  power  speetral  density  (i.e. 
noise  power  per  unit  bandwidth). 

Many  studies  have  shown  that  the  threshold  eontrast-energy  (E)  of 
a  target  is  linearly  related  to  the  power  speetral  density  (N)  of 
superimposed  broadband  noise  [1].  This  linear  relationship  is 
found  for  a  variety  of  targets  (e.g.,  letters  and  gratings),  although 
the  slope  of  threshold  energy  verses  power  speetral  density 
funetion  (or  E-N  funetion),  varies  with  the  type  of  target  used. 

In  this  study,  we  measured  energy  thresholds  on  both  a  speekle 
baekground  and  baekgrounds  eonsisting  of  randomly  modulated 
pixel  luminanee  (pixel  noise).  The  pixel  noise  data  were  then  used 
to  generate  E-N  plots  for  two  targets.  Further,  the  power  speetral 
density  of  the  laser  speekle  pattern  imaged  on  the  retina  was 
estimated  using  measured  speekle  eontrast  and  the  optieal 
parameters  of  the  human  eye.  We  found  that,  at  the  same  power 
speetral  density  levels,  energy  thresholds  obtained  for  gratings 
superimposed  on  pixel  noise  were  slightly  higher  than  those 
obtained  for  gratings  superimposed  on  speekle  noise.  That  is, 
speekle  noise  was  less  effeetive  than  pixel  noise  in  raising  energy 
thresholds. 


2.  Methods 

2.1  Observers 

Two  observers,  23  and  28  years  old,  partieipated  in  the  study. 
Both  observers  had  normal  or  eorreeted  to  normal  aeuity. 

2.2  Stimuli  and  Apparatus 

There  were  three  stimulus  eonditions  referred  to  here  as  speekle, 
uniform-luminanee,  and  pixel-noise.  In  the  speekle  eondition,  a 
grating  test-target  was  superimposed  on  a  baekground  speekle 
pattern.  In  the  uniform-luminanee  eondition,  a  grating  test-target 
was  superimposed  on  a  uniform,  baekground  field.  In  the  pixel- 
noise  eondition,  the  test-target  was  embedded  in  pixel-noise.  This 
image  was  then  superimposed  on  a  uniform  field.  In  all  eases  the 
grating  test-target  had  a  sinewave  luminanee  pattern  multiplied  by 
a  unit-height  Gaussian  aperture,  whose  standard  deviation  was 
about  0.48°.  The  width  of  the  target  was  defined  as  +/-  3  standard 
deviation  units  or  2.88  degrees.  The  spatial  frequeneies  of  the 
sinewaves  were  2.2,  4.4,  8.7,  and  17.4  eyeles/deg.  The  mean 
luminanee  of  the  test- targets  was  about  2.0  ed/m^.  The  grating 
test-targets  were  generated  in  Matlab  using  the  Psyehophysies 
Toolbox  extensions  [2,  3],  and  were  displayed  using  a  DEP 
projeetor. 

The  speekle  baekground  was  obtained  from  a  pre-produetion, 
seanning  laser  projeetor  (Evans  &  Sutherland),  and  the  uniform- 
luminanee  baekground  was  obtained  using  an  ECD  projeetor.  The 
size  of  the  relevant  baekgrounds  were  adjusted  so  that  they 
subtended  about  34°  horizontal  by  16°  vertieal  for  both  the  0.965 
and  1.93  m  viewing  distanees.  The  mean  luminanee  of  the 
baekgrounds  was  4.8  edW.  All  test  and  baekground  stimuli  were 
displayed  on  a  rear-projeetion  sereen  with  a  gain  of  1.0  and  a 
high-eontrast  tint  (Pro Sereen  Ine.). 

The  RMS  eontrasts  of  the  targets  and  baekgrounds  were  estimated 
from  images  eaptured  by  an  imaging  photometer  (Eumetrix  model 
500e).  The  eontrast  of  eaeh  image  pixel  was  ealeulated  by 
dividing  its  deviation  from  the  mean  image  luminanee  by  the 
mean  image  luminanee.  The  RMS  eontrast  was  eomputed  by 
taking  the  square  root  of  the  mean  squared  eontrast. 


Figure  1.  A  photograph  of  the  laser-speckle 
background  stimulus. 
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Figure  2.  A  photograph  of  the  pixel-noise  pattern 
that  was  added  to  the  test-stimulus  and  thus 
effectively  acted  as  a  background  image. 


2.3  Procedure 

Observers  first  adapted  for  8-10  min  to  the  ambient  illumination 
of  the  experimental  room.  Eaeh  experimental  trial  eonsisted  of 
two  temporal  intervals  demareated  by  tones.  Only  one  interval 
eontained  the  target,  and  the  observer  was  instrueted  to  push  a 
button  eorresponding  to  that  interval.  The  81%-eorreet  threshold 
was  estimated  using  the  Quest  proeedure  [4],  whieh  terminated 
after  30  trials.  The  four  target  spatial  frequeneies  were  tested  in  a 
random  order  for  eaeh  observer,  and  eaeh  observer  produeed  3-7 
threshold  estimates  for  eaeh  target/noise  eondition.  The  mean  and 
standard  error  of  the  mean  for  a  given  stimulus  eondition  were 
eomputed  using  the  eombined  data  of  both  observers.  Thus,  the 
eontribution  of  eaeh  observer’s  data  to  the  mean  was  weighted  by 
the  number  of  measurements  produeed  by  that  observer. 


3.  Results 

Figure  3  shows  RMS  eontrast  thresholds  as  a  funetion  of  target 
spatial  frequeney.  The  individual  eurves  represent  different  noise 
eonditions.  Speekle  inereases  eontrast  thresholds  relative  to  the 
uniform  (no  noise)  eondition  at  all  spatial  frequeneies,  eonsistent 
with  the  wide  spatial  frequeney  bandwidth  (small  speekle  size)  of 
speekle  noise. 

The  RMS  eontrast  of  the  speekle  noise  was  more  than  four  times 
that  of  the  pixel  noise — 0.29  and  0.07,  respeetively.  However, 
Figure  3  shows  that  speekle  noise  inereased  thresholds  less  than 
the  pixel  noise.  In  the  next  seetion  we  show  how  the  small  angular 
size  of  the  speekle  elements  reduees  its  effeetiveness  in  inereasing 
target  thresholds. 


3.1  Estimation  of  the  Power  Spectral  Density  of 
Pixel  noise 

If  independent  samples  of  zero-mean  Gaussian  luminanee  noise 
are  added  to  eaeh  square  region  (eheek)  on  a  display,  the  power 
speetral  density  is  defined  [5]  to  be  the  square  of  RMS  eontrast 
divided  by  the  two-sided  bandwidth  of  the  noise.  That  is. 


where  the  two-sided  noise  bandwidth  is  given  by: 
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Figure  3.  Contrast  threshold  as  a  function  of 
spatial  frequency. 


where  Ax  and  Ay  are  the  width  and  length,  respeetively,  of  the 
noise  elements  in  degrees  of  visual  angle. 

The  produet  AxAy  is  the  area.  A,  of  the  noise  elements,  and  it  is 
eustomary  to  eompute  N  as  follows  [2]: 

N=Cl,,-A  (de^),  (3) 

where  C  is  the  RMS  eontrast.  In  the  present  study,  the  RMS 
eontrast  of  the  pixel  noise,  measured  using  the  imaging 
photometer,  was  0.07.  At  the  1.93  m  viewing  distanee,  the  angular 
length  and  width  of  a  pixel  noise  element  was  57.40  •  10'^  degrees 
and  so  the  pixel  area  was  826.56  •  10'^  degrees^.  Using  these 
values  with  Equation  (3)  yields  a  power  speetral  density  of  4.05 
TO'^  degrees. 

For  the  0.965  m  viewing  distanee,  the  angular  width  and  length  of 
the  noise  element  double  whieh  inereases  the  angular  area  and 
power  speetral  density  by  a  faetor  of  four.  The  eontrast  and  area 
values  used  to  eompute  power  speetral  density  are  shown  in  Table 
1. 


Tablel  1.  Values  used  to  compute  PSD. 


Pixel  Noise 
1930  mm 
Viewing 
distance 

Pixel  Noise 
965  mm 
Viewing 
Distance 

Speckle 

4.31  mm 
Pupil 
diameter 

Contrast 

0.07 

0.07 

0.29 

Noise  element 
area  (lO'^deg^) 

826.56 

3306 

56.74 

PSD  (lO-^deg^) 

4.05 

16.2 

4.77 
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3.2  Estimation  of  the  Power  Spectral  Density  of 
Speckle 

The  width  of  a  speckle  element  on  the  retina  is  proportional  to  the 
pupil  diameter  of  the  observer.  Pupil  diameter  was  computed 
using  a  formula  developed  by  Barten  [6] 

p  =  5  -  3  tanhj  0.4 lo^Z  A /  40^ j  |  (4) 

which,  for  the  mean  luminance  L  (6.8  cdW  )  and  display  area  A 
(891  degrees^)  used  in  this  study,  results  in  a  pupil  diameter  of 
4.31  mm. 

The  average  width  of  a  speckle  element  on  the  retina  was 
computed  using  the  Airy  disk  formula: 

w  =  1.22aT  ,  (5) 

P 

where  X  is  wavelength,  /  is  the  distance  from  the  nodal  point  of 
the  eye  to  the  back  of  the  retina,  and  p  is  pupil  diameter  in  mm. 
Using  Equation  (5)  with  a  4.31  mm  pupil  diameter,  a  laser 
wavelength  of  0.532  pm  (the  green  channel  of  our  display),  and 
an /of  16.5  mm  yields  a  speckle  width  of  8.50  •  10'^  pm.  Because 

a  speckle  element  is  circular,  the  speckle  area  is  ;r(w/2)^  or  4.91 

pm^.  Finally,  the  angular  speckle  area  is  produced  by  multiplying 
the  speckle  area  in  pm^by  1.16  TO^degrees/pm.  [7] 

The  power  spectral  density  of  the  speckle  pattern  was  estimated 
using  Equation  (3).  The  contrast  and  area  values  used  to  compute 
this  power  spectral  density  are  shown  in  Table  1. 

3.3  E-N  Plots 

The  data  of  Figure  3  are  presented  in  Figure  4  as  threshold  energy 
verses  power  spectral  density  (E-N)  plots.  Energy  is  equal  to  the 
product  of  squared  RMS  contrast  and  stimulus  area  (6.15 
degrees^).  The  solid  lines  represent  linear  fits  to  the  two  pixel- 
noise  and  one  uniform-luminance  thresholds  (circles).  The 
speckle  thresholds  are  shown  by  the  square  symbols. 

The  speckle  thresholds  were  lower  than  thresholds  on  pixel  noise 
of  the  same  power  spectral  density.  For  the  2.18  cycles/deg 
grating  the  difference  was  marginally  significant  [t(10)  =  1.12,  p  < 
0.058].  For  the  4.36  cycle/degree  grating  the  difference  was 
highly  significant  [t  (10)  =  4.96, /?<  0.0003]. 

The  dashed  horizontal  lines  in  Figure  4  connect  the  speckle- 
threshold  data  points  to  the  function  fitted  to  the  pixel-noise  and 
uniform-luminance  threshold  data.  A  lines  dropped  from  this 
intersection  to  the  power  spectral  density  axis  provide  an  estimate 
of  the  equivalent  pixel  noise.  The  equivalent  pixel  noise 
represents  the  power  spectral  density  of  pixel  noise  that  produces 
the  same  threshold  as  speckle.  The  average  equivalent  pixel  noise 
was  1.97  •  10'^  degrees^  which  is  less  than  half  of  our  estimate  of 
the  power  spectral  density  of  speckle. 

4.  Discussion 

The  results  of  this  study  suggest  that,  at  the  same  power  spectral 
density,  speckle  may  be  slightly  less  efficient  than  pixel  noise  in 
increasing  target  thresholds.  One  possible  explanation  of  this 
finding  is  that  scattering  by  the  ocular  media  may  differentially 
affect  laser  and  non-coherent  radiation.  However,  before  we  can 


Figure  4.  E-N  plots. 


conclude  that  speckle  is  less  efficient  than  pixel  noise  in 
increasing  target  contrast,  improvements  to  the  method  used  to 
estimate  speckle  power  spectral  density  are  required 

Equation  (3),  which  was  used  to  estimate  power  spectral  density, 
is  only  accurate  for  a  Gaussian  luminance  distribution.  Shown  in 
Figure  5  are  the  luminance  distribution  of  the  pixel-noise  image 
(red  line)  and  a  Gaussian  distribution,  demonstrating  that  this 
condition  was  satisfied  for  the  pixel  noise. 

Shown  in  Figure  6  are  the  luminance  distribution  of  the  speckle 
image  (red  line)  and  a  Gaussian  distribution.  The  speckle 
luminance  distribution  is  clearly  not  Gaussian.  We  are  currently 
attempting  to  modify  the  equation  used  to  estimate  the  power 
spectral  density  of  speckle  in  order  to  account  for  the  non- 
Gaussian  luminance  distribution.  We  are  also  investigating  using 
discrete  Fourier  transform  techniques  to  estimate  the  power 
spectral  density  of  the  noise  patterns. 

It  is  clear  from  the  present  data  that  speckle  contrast  is  an 
incomplete  metric  for  characterizing  the  effect  of  speckle  on 
target  detection,  in  that  pupil  diameter  must  also  be  considered. 
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Pixel  Noise 


Figure  5.  Normalized  histograms  of  standardized 
luminance  calculated  as:  (Pixel  Luminance  - 
Mean  Luminance)/Standard  Deviation  of  the 
Luminance  Values. 


Speckle 


Figure  6.  Normalized  histograms  of  standardized 
luminance  calculated  as:  (Speckle-Element 
Luminance  -  Mean  Luminance)/Standard 
Deviation  of  the  Luminance  Values. 


Artigas  et  al.  [8]  used  artificial  pupils  and  also  found  that 
decreasing  the  diameter  of  the  pupil  increased  speckle-contrast 
thresholds. 

This  finding  that  pupil  diameter  must  be  considered  when 
assessing  target  detection  on  laser  speckle  has  practical 
implications  for  the  design  of  laser  displays.  As  the  luminance  of 
a  laser  display  increases,  the  pupil  diameter  of  observers  viewing 
the  display  will  decrease.  This  will  increase  the  power  spectral 
density  of  the  speckle,  making  targets  more  difficult  to  detect. 
Therefore,  methods  designed  to  reduce  the  speckle  contrast  [9] 
will  become  more  important  as  the  luminance  of  laser  displays  is 
increased. 
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